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Abstract 

While pure Yang-Mills theory feature the centre symmetry, this symmetry is 
explicitly broken by the presence of dynamical matter. We study the impact of the 
centre symmetry in such QCD-like theories. In the analytically solvable Schwinger 
model, centre transitions take place even under extreme conditions, temperature 
and/or density, and we show that they are key to the solution of the Silver-Blaze 
problem. We then develop an effective SU(3) quark model which conhnes quarks by 
virtue of centre sector transitions. The phase diagram by confinement is obtained 
as a function of the temperature and the chemical potential. We show that at low 
temperatures and intermediate values for the chemical potential the centre dressed 
quarks undergo condensation due to Bose like statistics. This is the Fermi Einstein 
condensation. To corroborate the existence of centre sector transitions in gauge 
theories with matter, we study (at vanishing chemical potential) the interface tension 
in the three-dimensional gauge theory with Ising matter, the distribution of the 
Polyakov line in the four-dimensional SU(2)-Higgs model and devise a new type of 
order parameter which is designed to detect centre sector transitions. Our analytical 
and numerical findings lead us to conjecture a new state of cold, but dense matter 
in the hadronic phase for which Fermi Einstein condensation is realised. 
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1 Introduction: 



A great deal of efforts, theoretical and experimental, are devoted to explore the properties 
of hadronic matter under extreme conditions, temperature and/or baryonic densities. The 
general belief is that matter is organised in a quark-gluon plasma phase in this regime. At 
least for small densities, this belief is corroborated by lattice gauge calculations [lH3] and 
collision experiments such as undertaken at RHIC [HIS]. 

Central to understand matter under extreme conditions is the understanding of colour 
confinement since its realisation is, almost by definition, the most prominent difference 
between the so-called hadronic and the quark-gluon plasma phase. Early on in the eighties, 
it was pointed out that the centre of the gauge group plays a major role for confinement: 
Realisation of the centre symmetry of the gauge sector goes in line with confinement 
while its spontaneous breakdown (e.g. at high temperatures) signals colour liberation [6l[7]. 
Roughly at the same time, it was proposed that certain degrees of freedom of Yang-Mills 
theories, such as monopoles or vortices, are responsible for the (dis-)order of the centre 
symmetry (see [8] for a review). It took until the late nineties to isolate these degrees of 
freedom in lattice gauge theories in a physical, i.e., regulator independent way PHU]. It 
was pointed out in [TT] that the mere existence of centre vortices as physical degrees of 
freedom demand a fine tuning between vortex entropy and energy. Quite recently [12], 
smooth stream-line configurations which bear confinement have been found using lattice 
gauge simulations: it was pointed out that if centre vortices are images of these stream- 
line configurations in a certain gauge, it would naturally account for the intrinsic fine- 
tuning. Most importantly, the centre vortex picture offers an understanding of confinement 
on the basis of weakly interacting degrees of freedom, which explain high temperature 
deconfinement and the centre disorder of spatial correlations at the same time [T3HT5] . 

Although a quite detailed picture of confinement and deconfinement at finite temperatures 
has emerged over the last decade, very little is known about cold but dense matter from 
first principles. Because of the notorious "sign" -problem, Monte-Carlo simulations of Yang- 
Mills theories with dense quark matter are only feasible for SU(2) colour [16] for which the 
sign problem is absenlll] or for rather small densities (see [IB] for a recent review). First 
insights into the properties of matter which feature in cold but dense QCD might be gained 
from exact solutions of models which mimic certain aspects of QCD. As an example, we 
mention the Gross- Neveu model which features the baryon crystal as an hitherto unknown 
state of matter [T9142T] . 

Among the very few gauge theories which admit an exact solution even for the case of 
dense quark matter is 2d Quantumelectrodynamics, the so-called Schwinger model [22] . 
The model with massless fermions was exactly solved in Hamiltonian formalism on the line 
in [231125] and on in [261127]: chiral symmetry is spontaneously broken and only states 
with a vanishing net baryon number appear in the spectrum. The model on the torus has 

Yang-Mills theory with dynamical fermions has no sign-problem as well and simulations at finite 
temperature and density are under way [TT] . 
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been studied in [28] and in particular the temperature dependence of the chiral conden- 
sate, Wilson loop correlators and Polyakov line correlators have been determined [29143T] . 
In [32l[33], non- vanishing values of the fermion chemical potential have been firstly con- 
sidered. It was found, most importantly, that the full non-perturbative partition function 
is independent of the chemical potential. Since then, many more quantities have been ob- 
tained in the massless Schwinger model ranging from the temperature dependence of the 
correlators of hadronic currents, spectral functions to the screening mass [3l] . This makes 
the Schwinger mode the ideal testbed to test new ideas, and we will make extensive use of 
this below. 

Besides of exactly solvable models, understanding mechanisms, which have been revealed 
in model studies or by non-perturbative approximations, is an invaluable tool since they 
might extend their applicability to the theory of interest, in particular, cold and dense 
QCD matter. Recent examples are the proposal of the quarkyonic phase the existence 
of which has been motivated by the large Nc considerations {Nc being the number of 
colours) [351 136]. Another recent example is the chiral magnetic effect which describes 
an induced electromagnetic current alongside an external magnetic field made possible by 
topological charge transitions in the quark gluon plasma phase [371138] . 

Of particular importance for this publication is the recent observation that quarks effec- 
tively comply to periodic boundaries conditions if exposed to a particular non-trivial centre 
sector in SU (Nc) Yang-Mills theory for A^^ even. In the dense hadronic (confining) phase, 
this opens the possibility that centre dressed quarks undergo condensation reminiscent of 
Bose- Einstein condensation [39]. In analogy, this has been called Fermi Einstein conden- 
sation (FEC) [10]. We also point out that the sensitivity of the quark spectrum to the 
boundary conditions of the quarks has been studied in [H]. 

In this paper, we will further study the phenomenological impact of the FEC effect. We 
will use the Schwinger model which will allow to study FEC on almost purely analytical 
grounds. Since in previous studies we considered an even number of colours, we here extend 
the considerations to the more relevant case of an SU {N^ = 3) gauge group using a quark 
model. We are going to show that transitions between the centre sectors are sufficient 
to confine quarks in this model. The quark model itself determines the so-called centre 
sector weights which serve as an order parameter for confinement. Using these weights, 
we will be able to calculate the phase diagram of the model as a function of the chemical 
potential and the temperature. For Nc = 3, we find that FEC does occur under cold and 
dense conditions at the presence of pressure. Essential for FEC is that the gluonic states 
manoeuvre through the centre sectors. Since non-trivial centre-sectors are suppressed by 
the presence of dynamical matter, we accumulate evidence in the remainder of the paper 
that centre sector transitions do occur in the so-called hadronic phase. To this aim, we will 
study the gauge theory with Ising matter and the SU(2) Higgs theory. We develop a 
novel order parameter which is sensitive to centre sector transitions and provide numerical 
evidence that these transitions take place until centre symmetry is spontaneously broken 
at high temperatures. 
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The paper is organised as follows: In section [21 we extend the considerations of |12] - H^ and 
show that the Yang-Mills "empty" vacuum possess gauge in-equivalent flat directions which 
collapse to the so-called centre-sectors once quantum fluctuations are included. Section |3] 
addresses FEC in the context of the Schwinger model. In section HI we develop the SU{3) 
quark model with confinement by virtue of the interaction of the quarks with the centre 
sector background fields. Confinement is established by studying the model's thermal 
excitations, and the phase diagram from confinement is obtained. In section |5l we start to 
investigate centre-sector transition by means of numerical simulations. For this purpose, 
we study the gauge theory with Ising matter, which has not yet been studied before 
to our knowledge. In section El we extend the simulations and consider the SU(2) Higgs 
theory. The order parameter for centre sector transitions is developed and our numerical 
findings for this order parameter are presented. Conclusions are left to the final section. 



2 Yang-Mills moduli space and centre sectors 
2.1 The empty vacuum on a torus 

A configuration with minimal Euclidean action, often called empty vacuum, is often the 
starting point of perturbation theory. In Abelian theories, it is naturally defined as a state 
for which the field strength at any point of space-time vanishes. In non- Abelian theories, 
such as SU{Nc) Yang-Mills theories, a more stringent definition is in order: there, the 
empty vacuum is a state for which the holonomy calculated along any contractible loop C 
yields the unit element of the gauge group: 

P exp |i J A^{x) dx'^^ = 1 (empty vacuum condition), (1) 

where A^(x) is the gauge potential and P denotes path-ordering. Throughout this paper, 
we will consider a 4-torus as space-time manifold. The gauge potentials Afj,{x) and gauge 
transformations ^l{x) satisfy periodic boundary conditionqj, matter fields such as a scalar 
Higgs field (f){x) or quarks q{x) are subjected to periodic and anti-periodic boundary condi- 
tions, respectively. We will adopt a lattice regularisation with the lattice spacing a acting 
as an UV regulator. Thereby, the gauge degrees of freedom are represented by the links 
U^{x) G SU{Nc). Gauge transformations act as usual: 

U^{x) = Q{x)U,ix)Q^{x + fi) . (2) 

In lattice discretisation, the vacuum condition ([TD becomes 

Y[Ui = 1, ^=(x,/i), (3) 

eec 

where C is a closed and contractible path on the lattice. The smallest contractible loops 
on the lattice are the loops surrounding the elementary plaquettes of the lattice. For the 

^In the continuum theory this implies a vanishing instanton number. 
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Figure 1: Absence of Polyakov line correlations in the "empty" vacuum (left). Step 1 of 
the complete gauge fixing (right). 



later choice, the vacuum condition ([3]) implies that the non-Abelian field strength vanishes 
for the empty vacuum. 

Already in the late nineties, it has been discovered that there are smoothly connected 
configurations which are gauge in- equivalent and which all satisfy the vacuum condition ([T]) 
(or ([3])) [I2HII]. The set of these configurations defines the Yang-Mills moduli space. While 
Keurentjes et al. discuss 3-dimensional Yang-Mills theory, the arguments were extended to 
4-dimension in [U]. To make the paper self-contained, we will give a full account of the 
Yang-Mills moduli space using lattice regularisation. 

Of particular importance for the study of confinement in Yang-Mills theories on the torus 
is the Polyakov line: 

V{x) = l[Uo{t,x) . (4) 

t 

Let us firstly study their correlations in the empty vacuum. For this purpose, we consider 
the maximal (contractible) loop in figure [H left panel (note that the product V of spatial 
links is the same at the lower and upper time-slices due to periodic boundary conditions). 
The vacuum condition implies 

V{x)VV\y)V^ = 1 tiVix) = tiV{y). (5) 

Hence, the trace of the Polyakov line is necessarily constant for any choice of an "empty" 
vacuum. Given that the static quark anti-quark potential V{r) can be extracted from the 
Polyakov line correlator: 

{P{x) P{y)) oc exp {-V^(r)/T} , P{x) = ^ tr V{x), (6) 

where r = \x — y \ and T is the temperature, the finding ^ implies that a single "empty" 
state cannot sustain quark confinement. 



6 



1 














1 














1 














1 














1 














1 


u 












p ' 








Pk 


1 Pi, 




p 




p 












u 


1 


1 


1 


1 


1 




u 


1 


1 


1 


1 


1 



Figure 2: Gauge fixing step 2 (right) and (almost) completely gauge fixed configuration 
(left). 

Genuinely different vacua states are obtained by identifying all states which are related 
by gauge transformations. In order to calculate those, a complete unique gauge fixing is 
in order. In a first step, we adopt a Polyakov type of gauge fixing: all time like links 
Uo{t,x) for t > 1 are gauged to the unit element (see figure [U left panel). The only 
non-trivial time-like links are those at t = 1, which thus equal the Polyakov line. The 
gauge transformation Qi, Q2, etc, are not used for this first gauge fixing step. These gauge 
transformations are now employed to gauge almost all spatial links at time-slice t = 1 to 
the unit element (see figure |2], right panel). The empty vacuum condition then implies 
that all spatial links without arrows in figure |2] are transformed to the unit element. The 
spatial links U which are far most to the left then equal the spatial Polyakov line in this 
direction. We then again use the vacuum condition for the contractible loop in figure [2l 
left panel and conclude that 

Vixk) 1 V\x,) 1 = 1 ^ V{xk) = V{x,) = V . 

In this gauge and for an empty vacuum state, not only the trace of the Polyakov is constant 
but also the full Polyakov line including off-diagonal parts. After the above gauge fixing, 
only gauge transformation Vli remains unfixed and acts as a residual (7/060/ colour symmetry. 
In a final step, we consider the plaquette in the bottom left corner of the lattice in figure [21 
right panel. Generically, the eigenvalues of the spatial and time-like Polyakov lines are 
different. In order to satisfy the vacuum condition for this case, the matrices V and U are 
drawn from the Cartan subgroup of SU {N^ implying: 

VUV^U^ = 1 [p,U]=0 . 

Let the spatial Polyakov line U now be an element of the Cartan subgroup with all eigen- 
values different. The Polyakov line V must then be an element of the Cartan subgroup in 
order to describe an empty vacuum: 
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For two choices of V with at least one eigenvalue being different, we obtain two gauge 
in-equivalent empty vacua states. Note that the Polyakov line homogeneously transforms 
under gauge transformations, 

implying that its eigenvalues are gauge invariant. Each U{1) is spanned by a compact 
angle variable On, n = 1 . . . — 1, which span the space of global minima of the classical 
Yang-Mills action - the moduli space. For each direction, the corresponding Polyakov line 
can be chosen from the [f/(l)]^=~^ subgroup. Hence, the moduli space is at least spanned 
by the 

group manifold. Note that if, for a particular direction, the eigenvalues of the Polyakov 
line are / fold degenerate, the Polyakov lines of the other direction can be chosen from a 
SU (/) subgroup and still satisfy the empty vacuum condition. Hence, the moduli space is 
slightly larger than the space spanned by the U{1) groups only. 

We finally make two comments: 

(i) We point out that the trace of the Polyakov line being different is a sufficient, but not a 
necessary condition for two states being at different points in moduli space: assume that 
two Polyakov lines possess different eigenvalues but the same sum of all eigenvalues. They 
would belong to different points in moduli space, but their trace would be equal. 

(ii) A perturbative treatment should involve a summation over all states with minimal 
(global) action. This implies an integration over the moduli space. Note that the standard 
perturbation theory merely chooses one state (i.e., the state of vanishing gauge potential or 
the state of all unit links in lattice formulation) of the moduli space. The integration over 
the moduli space would most likely remove the colour states from the theory thus inducing 
confinement. It would not, however, provide a confinement scale. Perturbation theory 
with an integration over the moduli space as well as the phenomenological implications of 
the moduli space integrations in effective quark theories will be explored elsewhere. 



2.2 Yang- Mills quantum vacuum 

The flat directions of the "empty" vacuum, discussed in detail in the previous subsection, 
are lifted if quantum fluctuations are considered. The symmetry of the quantum effective 
action collapses to the discrete centre symmetry. Introducing the centre elements of SU{Nc) 
by ^ 

Zm = expji-^mj, m = 1 . . . N^, (7) 

the centre transformed lattice configuration {f/^} is obtained by multiplying all time-like 
links at a given time slice to by Zm- 

f/o(to,a;) = ZmUo{to,x) , for Va; , (8) 
U'^it,x) = U^it,x), else. (9) 
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If the lattice action in pure Yang-Mills theory consists of a collection of contractible loops 
(such as the Wilson action which is constructed from the smallest of such loops - the 
plaquette), the action is invariant under the above centre transformation. On the other 
hand, the Polyakov line (jlj) transform homogeneously: 

V[U%x) = z^V[U]{x) . (10) 

Any ergodic Monte-Carlo simulation on a finite lattice necessarily averages over the centre 
copies in a democratic way implying that 

{P[U]) = (P[t/1) = {P[U]) {P[U])=0, 

which is also true in the high temperature phase of Yang-Mills theory. Showing that centre 
sector transitions imply confinement hence request more subtle arguments involving centre 
invariant expectation values. 

To this aim, let us directly consider the static quark antiquark potential V{r) as inferred 
from the Polyakov line correlation function ([6]). Instead of the Polyakov line expectation 
value, we consider its spatial average in relation to a reference P{0) on the lattice. We 
observe 

(^P(O) ^P(a;)^ oc ^e-^('^)/^ = finite lim V{r) ^ oo , 

X X 

and, hence, confinement. If on the other hand the potential approaches a finite value for 
r — )■ oo, the correlator necessarily behaves as 

lim (P(0) P(a;)) = finite or (p{Q)^P{x)^ ^ oo . (11) 

X 

In view of f lTOj) . the latter equation could mean that centre sector disorder does not occur 
at length scales set by the lattice size. In the next subsection, we will show that the latter 
condition is in line with what is usually referred to as the spontaneous breakdown of centre 
symmetry. This breakdown occurs at high temperatures leaving us with the so-called 
quark-gluon plasma phase. 

2.3 Spontaneous breaking of centre symmetry 

One way to reveal the spontaneous breaking of centre symmetry using ergodic Monte-Carlo 
simulations is to add a centre symmetry breaking source term to the action and controlling 
the strength of this term by a parameter, let us say, /i, which is reminiscent of a magnetic 
field in an Ising type setting. To be explicit, we will study SU(2) lattice gauge theory using 
a Euclidean space-time lattice represented by a A^^ x Nt grid with lattice spacing a. Using 
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SU(2), N x4, 1020configs 
1 I I I I I I I I I I I I I r 




Figure 3: Sketch of the response of the Polyakov hne expectation value to external centre 
symmetry breaking (left). The response function of the Polyakov line to an external field 
as function of /3wii for several spatial volumes (right). 

the Wilson action, the partition function is given 

Zym = j T>U exp{S'wil + 5'break} , (12) 

^Wii = ^ E tT{U,{x)U,{x + i,)Ul{x + u)Ul{x)] , 

5'brcak = ^ X] -P(^) ' 

X 

where we have used the Polyakov line (jl]) to test the response of the theory to explicit 
centre symmetry breaking. 

In the confinement phase, the Polyakov line expectation value (P) linearly responds to the 
presence of an external centre symmetry breaking. This is illustrated in figure [3l left panel, 
curve (1). Let us now consider the quark gluon plasma, i.e., the deconfinement phase. For 
finite lattice volume, the Polyakov line expectation value necessarily vanishes at h = 
and is otherwise a smooth function of h. This is illustrated by curve (2) in figure [31 left 
panel. Only in the infinite volume limit, the (P) becomes a discontinuous function which 
approaches a finite value in the limit /i — >■ (see curve (3) in figure [3l left panel). 

How can we anticipate the spontaneous breaking of the centre symmetry by means of a 
finite volume ergodic Monte-Carlo simulation? 
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A straightforward way is to study the gradient of {P){h) aX h = 0. This gradient should 
rise beyond bounds with increasing volume. Hence, we find 



d{P){h) 



Imi (p(0)5^P(:r)) - (p(0))(5^P(:r)) 



dh 



h=0 



X X 



(p(0)5^P(a:))[^^ ^ oo. 



(13) 



X 



This agrees with the idea of ceasing centre sector transitions in the deconfinement phase 
(see (fTT]) ). Figure [3l right panel, shows the response function ( IT3|) for a A^^ x 4 lattice as 
a function of /Swu for several values of N. For this lattice geometry, the deconfinement 
phase is attained for /3 > P^cc ~ 2.3. We observe that the response function is basically 
independent of the spatial size in the confinement phase for /3 < /3dcc while it roughly scales 
with the volume in the gluonic plasma phase. 

3 Fermi-Einstein condensation in the Schwinger model 



It was already observed in [39] that centre sector transitions have far reaching phenomeno- 
logical consequences for gauge theories with dynamical fermions at finite densities (such 
as QCD): for an even number of colours, the centre transformed background field can be 
viewed as imposing periodic boundary conditions for the matter fields. This opens the possi- 
bility that the centre dressed fermions undergo condensation similar to that which is known 
as Bose-Einstein condensation although the fermions are described by anti-commuting 
Grassmann fields. In analogy, this new phenomenon has been called Fermi- Einstein con- 
densation [10]. For this scenario, centre transitions are essential hence relegating Fermi- 
Einstein condensation to the confining phase at intermediate values of the chemical po- 
tential. Consequently, Fermi-Einstein condensation evades the spin-statistic theorem since 
the fermion fields cannot be viewed as asymptotic states. 

Here and in the next section, we further explore Fermi-Einstein condensation: (i) the 
U(l) gauge theory with fermions at finite densities, i.e., the Schwinger model, bears all 
prerequisites for Fermi-Einstein condensation. Because it can be solved analytically, it is 
also the ideal testbed to work out the consequences of the condensation, (ii) The rise of 
Fermi-Einstein condensation is QCD-like theories with an odd numbers of colours is less 
clear cut. We therefore investigate a QCD inspired SU(3) effective quark model to explore 
its possible existence and phenomenological consequences. 

3.1 Schwinger-model essentials 

Here we illustrate the role of the centre with the finite temperature Schwinger model in 
a spatial box of length L. The fermion field is anti-periodic and the gauge potential 
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periodic in Euclidean time with period /3 = 1/T. Only configurations with vanishing 
instanton number contribute to thermodynamic potentials such that we may assume that 
the fields are periodic in the spatial direction. Thus we consider the Schwinger model on 
the Euclidean torus [0,/3] x [0,L] with volume V = f3L. Only at the end of the calculation 
do we let the spatial extension L tend to infinity. 
To continue we decompose the gauge potential as |29] 

Ao = -jho + doX - di(j), = —h + diX + (14) 

with constant toron fields and hi. In the zero-instant on sector the mapping {/i^, 0, A} — )■ 
Afj^ is one-to-one if we assume that the periodic functions A and (f) integrate to zero. Adding 
an integer either to Hq or to hi is equivalent to performing a gauge transformation with 
winding and thus we must further assume that < < 1. The Polyakov-line is given by 



P{xi) = exp |i J AodxQ^ = exp{27ri /iq} exp |— i9i J (j){x) dxc 



(15) 



Using the [0, 1] periodicity of the partition function in Hq, the shift 

ho ^ ho + a (16) 

transforms the Polyakov line expectation value as 

(P) ^ exp{27ri a} (P) , aE [0, 1] (17) 

and is identified as U{1) centre transformation of the gauge fields. Hence, integration over 
the toron field ho implies an average over the centre sectors of the theory. 

The Jacobian of the transformation (fT4|) can be calculated by expanding all fields in Fourier 
modes and this way one finds [29] 

VA^ = {27iYdet'{~A)V(j)VXdhodhi, (18) 

where the primed determinant is the product of all positive eigenvalues of the operator —A 
on the torus. The functional integral over the gauge functions A cancels out in expectation 
values for gauge invariant objects and the 0-dependence of the fermionic determinant can be 
calculated explictly by integrating the chiral anomaly |15] . The resulting functional integral 
over the periodic function (which must integrate to zero) is Gaussian and leads to the 
following expression for the grand canonical partition function Z(/3,L,/i) = tr ^Q-P^+^iNy^ 



Z{/3,L,f,) = dett-ifL) Z"^^" ^'^'^^ ^^^^ 

z{ho; /3,L,fi) = [ dhi det{i$h,^) , (20) 
Jo 
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where = el\fv: is the induced photon mass, and zih^] l3,L,fi) is the centre sector 
partition function specified by a particular value of Hq. The Dirac operator is given by 

= i$- (^jho- i^) 7° - (^/^i) 7^ . (21) 

The fermionic determinant det(i^/j^^) is an elliptic function of the constant toron fields 
ho, hi and depends on the chemical potential. We use anti-periodic boundary conditions 
in time direction and periodic ones in spatial direction. 



3.2 Centre-sector partition function 

Let us further study the grand-canonical partition function z{ho; P, L, fi) for a given centre 
sector. This study emulates the spontaneous breaking of centre symmetry where the theory 
is subjected to a "frozen" value of ho. We stress that such the spontaneous breakdown 
cannot occur in 2 dimensions due to severe infra-red divergences induced by the would-be 
Goldstone bosons. We will however find that the suppression of centre sector transitions 
will lead to the wrong physics, and we believe that this is a generic truth valid also in 
higher dimensions. 

To recover the familiar Fermi sphere physics, it is sufficient to consider large spatial volumes 
for the moment. In this limit, we will be able to analytically perform the integration over 
the spatial component hi. To start with, we find: 

det(i4,^) = e-PLfif.,p,LM) ^-PEc^LM) (22) 

where Ecas{L,hi) is the independent Casimir energy (details of the calculation are left 
to appendix R|) : 

^cas = -— + — (l-2hi)\ (23) 
The /X and temperature dependence is encoded in the so-called 'free energy density' 

/ = -^^ln[(l + e2-'^«e-^(^"+'^)) (i + e-^-^Og-zsCi^n-/.)) j ^ (24) 

where L En = 27T\n — hi\. For — )■ oo the free energy density does not depend on the 
constant gauge field hi. In a theory with fixed centre sector, i.e., with a frozen value ho, 
the baryon number density would be given by 

= l/n\ = _^ 1 f 1 1 1 ...^ 

L^"^' dfi TiJo ^\e2-''oe/3(p-M) + i e-2.i''oe/5(p+M) + 1 / • ^^^^ 

For the trivial centre sector ho = 0, we obtain the standard result of a free Fermi gas 
forming a Fermi sphere at finite values for /i. On the other hand, for the sector ho = 1/2, 
we obtain a scenario as if the fermions were to acquire Bose-statistics: 

1 poo 

/(/X, ho = 1/2) "^^-—j^ dp[ log (1 - e-^^^-'^)) + log (1 - e-^(P+^)) } , 
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Note that the latter free energy logarithmically diverges when p — )■ This is the usual 
instability due to condensation and corresponds to Fermi-Einstein condensation in the 
present context. Note that this singularity is integrable upon the integration over the 
centre sectors Hq. 

In the remainder of this subsection, we will show that the assumption of a fixed centre sector 
will lead to the wrong physics in the case of the Schwinger model. The key observation is 
that the only excitations of the model are fermion anti-fermion bound states. The theory 
lacks any physical states which would couple to the chemical potential and would give rise 
to non-vanishing fermion density. The baryon density ought to vanish even in the case of 
non- vanishing values of the chemical potential. An inspection of fl25|) shows that this is not 
the case if we consider a fixed value of Hq only. The momentum integral can be performed 
leaving us with (T = 



n pb = fi + T 



In (1 + e"'"''^" e^/^) - In (l + e^"^'^" e^/^)] , ho ^1/2. (26) 



At /i = 0, the contributions from fermions and anti-fermions to the baryon density cancels 
since there is no fermionic mass gap in the theory. For ho = 1/2, the theory produces a 
condensed states, and the result ( 1261) cannot be extended to cover this case. For the trivial 
centre sector ho = 0, the baryon density linearly rises with the chemical potential. This 
is the well known Fermi-sphere behaviour for the 2-dimensional case. For any other value 
of ho, the baryon density is complex which renders the findings difficult to interpret in 
physical terms. 

The dependence of the baryon density on the chemical potential for a fixed centre sector 
is clearly spurious since the spectrum of the theory is free from states which carry baryon 
charge. A similar spurious dependence has been encounter in other theories as well which 
have been treated in an approximative way, and has been called the silver blaze problem [16]. 



3.3 Centre sector average 

Considering individually the contributions of the centre sectors (specified by ho) the free 
energy density has led to a singularity for ho = 1/2 associated with the condensation of 
fermions by virtue of periodic boundary conditions in this specific sector. In the Schwinger 
model, we are now in the comfortable situation that we can explore the physics of this 
Fermi-Einstein condensation since the centre sector average can be performed analytically. 
To this aim, we point out that the average of the fermion determinant over the spatial 
toron field can be written as (see appendix IA.2I for details): 

dh, det(i4,,) = E e-^-^^^+^-'^^ , ^ = ho + ffi, (27) 

Jo V2r|r/(ir)|2 ^ 2n 

where ?7(ir) is the Dedekind eta-function (see f llOip ). The centre sector average, i.e., the 
integration over ho, can be easily performed and restricts the Matsubara momentum sum 



14 



in ( p7|) to the trivial momentum p = 0: 



dhodhi det(i4,^) = ^f^j^^^ ■ (28) 

The Roberghe- Weiss symmetry for a model with global U{1) centre-symmetry implies that 
the partition function can not depend on the imaginary part of fi. Analyticity suggests 
that this might also hold for real values of the chemical potential. Indeed for the present 
case, we find that any dependence on disappears after the integration over ho. 

From [17] we take the result 

det'^/'(-A) = /3|r/(ir)|2 (29) 



VI 

dhodhi det(i(^^,^) = y ^^^/i/2/_^n - (^0) 



to write 



Inserting the latter finding into the grand canonical partition function f[T9|l . the square 
root of the determinant cancels, and we finally obtain 

Z{/3,L,fi) = J\ ^ ^ (31) 



2 ,/det'f-A 



2^ 



The key observation is that the /i dependence has disappeared implying that the centre 
sector average has solved the silver blaze problem: the baryon density vanishes independent 
of the value for the chemical potential as it should for the Schwinger model. 

It is instructive for QCD model building to investigate how the centre sector average 
eliminates the //-dependence. To this aim, we reconsider from the partition function f[T9|) 
the factor 

[' dho det(i4^) = e-^^ca.(L) f\f^^ -Q (1 + e2.i/^Oe-/3(^n+^')) (^i ^ Q-2.iho^-(3{E„-|.)^J ^ 

(32) 

where L En = 27r|n — hi\. Expanding the product yields terms such as 

exp{27rn i/io} exp{—(3nfi} . 

These terms vanish upon the integration over ho unless we have n = 0. Hence, the integral 
f p2|) does not depend on /i. This finding has a direct physical interpretation: n = only 
occurs if as many fermions as anti-fermions contribute to the product. Hence, the ho 
integration eliminates all states which carry a net baryon charge thus solving the silver 
blaze problem. 
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hO hO 



Figure 4: The real part of the differential contribution z in (120 p to the fermion determinant 
as a function of and for /3/i = 4 (left). The real part of 2 as a function of /iq and 
for = 1 (right). 

3.4 Volume studies 

The question whether a particular centre sector is singled out or whether a more democratic 
average over the centre sectors is in order should be answered by the theory itself. For 
this investigation, we will here consider the large volume limits L — )■ oo and /3 — )■ oo and 
combinations of it. 

In ( IT9l) . we decomposed the grand canonical partition function Z into a part from dynamical 
fields and a part z{}iq\ [5/L,[5^) featuring the contributions from the toron fields. From 
the studies in the previous subsection, we already know that Zfcr is actually independent 
of /3 /i. For a given value of the fugacity yu/T = we have calculated 2; as a function 
of /lo and P/ L. The colour coded result is shown in figure HI left panel. The important 
observation is that the largest contribution arises from the trivial centre sector around 
/io = at least for (3fi > 1 (note that the z is periodic, i.e., z{ho) = ziho) + n, n integer). 
The maximum at ho = even diverges in the limit (3/L — )■ 00. 

Here is a lesson to learn if it comes to Monte-Carlo simulations. Assume that we were to 
estimate the ho integral using a Metropolis Monte-Carlo method. First of all we note that 
for Hq ^ {0.5, 1} the determinant is complex. For small enough let us say < 1, 
we could use the real part of the determinant as a reweighting factor for the simulation. 
From figure IH left panel, it is already clear that the method samples the region around 
/lo = to a large extent. For /9/L ^ 1, other contributions from larger values of 
would be extremely rare. If we now stick to fixed aspect ratio we will find that for 
large values of the real part of the determinant develops negative parts invalidating 
the Monte-Carlo approach altogether (see figure HI right panel). The keypoint, however, is 
that the integral of z over /iq yields a constant only dependent on the aspect ratio (see the 
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Figure 5: The Polyakov line expectation value for vanishing chemical potential as a function 
of the temperature T and spatial size L (left). (P) for a fixed aspect ratio /3/L = 1 as a 
function of fi and T (right). 

previous subsections). This means that integrating along the horizontal axis in figure HI 
right panel, must produce the same value for every choice of (3fj,. Obviously, sampling 
the whole range of Hq values in the Monte-Carlo simulation is of key importance although 
difficult to achieve for P/L ^ 1, i.e., in the thermodynamic limit L — )■ cxd with T ^ fixed. 
Our findings suggest that the infinite volume zero-temperature limit of the QFTs (which 
we are going to discuss below) is delicate: we suggest that these limits should be taken at 
fixed aspect ratio t = (3/L and with one of the variables, L or /3 tending to infinity. In this 
context, the value r = 1 is of particular relevance for lattice gauge theories, since in this 
limit the rotational symmetry is recovered in the scaling limit of vanishing lattice spacing. 

Although the centre transitions never cease in the Schwinger model, the trivial centre sector 
contributes an overwhelming part to the expectation values. This is also clear from the 
Polyakov line expectation value. This expectation value can be calculated in closed from 
for arbitrary lengths /3, L and at presence of chemical potential fi (details of the calculation 
can be found in appendix IA.3P : 

(P) = e'^^^ exp (-^coth ^^^^ , (33) 

where = e/-\/vr is the dynamical generated photon mass which sets the scale of the 
theory. 

Let us firstly discuss the case of vanishing chemical potential. Our findings for this case 
are summarised in figure [5], left panel. Note that the Polyakov line expectation value (P) 
always vanishes for any fixed spatial size L in the zero temperature limit T — 0. On 
the other hand for a fixed temperature, (P) increases with increasing system size L, and 
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potentially approaches quite large value in the infinite volume limit: 



lim (P) = 

L— >oo 




We will see below that this behaviour will be mirrored by the SU(2) Higgs theory. 

Secondly, we consider the case of a non-vanishing chemical potential. For a fixed aspect 
ratio = 1, figure [5l right panel, shows the Polyakov line expectation value as a function 
of n and T. We point out that the fermion determinant is complex implying that the 
expectation value of a phase, in particular the Polyakov line, is not necessarily bounded 
anymore. These imaginary parts allow for the cancellations which are essential to solve 
the Silver Blaze problem and, at the same time, bear the potential for ^ 1. In the 
infinite volume limit, we now find 

hmjP) = exp|-^^| , 

which is bigger than one if the chemical potential exceeds the photon mass. 



4 Fermi-Einstein condensation in a QCD quark model 
with confinement 

4.1 Model building 

In order to trace out the phenomenology of centre sector transitions and Fermi-Einstein 
condensation in a more QCD-type setting, we are here going to investigate an effective 
SU(A'^c = 3) quark theory in four dimensions. The so-called constituent quarks q{x) satisfy 
the usual anti-periodic boundary conditions in Euclidean time direction (and periodic ones 
in the spatial directions) 

q{x + (3eo) = (-l)g(x), q{x + /3ek) = q{x), k = 1,2,3. (34) 

The quarks possesses the constituent quark mass m and are assumed to only interact with 
an homogeneous temporal gauge field specifying the centre sector: 

aJ"^ = 27mTH, 1 <n<N^, (35) 

where the generator H is from the Cartan algebra, i.e., H = diag(l, . . . , 1, 1 —Nc)/Nc. The 
Polyakov line is in the centre of the group and the centre sector is labeled by n since the 
trace of the Polyakov line is given by 

P._l,,exp{./^.„4"'}^.„. .„^e.p{?^.}. (36) 
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The crucial observation [39] is that by means of a Roberghe- Weiss transformation quarks 

^0 



which are subjected to the background field aI^^ can be considered as quarks manoeuvring 



in the trivial background A^^"^ but with changed boundary conditions: 

qix + /3eo) = -g(x) ^ g(x + /3eo) = -^„g(x) . (37) 

For even N^, there is the element zi\f^/2 = ~1) and it is this centre sector which is mapped 
to the trivial sector but with quarks now obeying periodic boundary conditions. Since it is 
essentially the boundary conditions which dictate the behaviour of the thermodynamical 
potentials, centre-dressed quarks follow Bose statistic and undergo condensation if the 
chemical potential approaches the fermionic mass gap [39] • Apparently, this scenario relies 
on the fact that (—1) is an element of the centre of the gauge group. This is only the 
case if the number of colours Nc is even. The more interesting, i.e., the more QCD type, 
case is A'^c = 3 and evades the line of arguments. In this subsection, we will explore the 
phenomenology of the centre sector transitions in the realm of an effective quark model for 
iVc = 3. 

The partition function of our model is given by 



/ DgDg exp|g(i^+ (aJ"^ + i/i)7° + im)g| = ^e 

n=l n=l 



p(n) 



(3^ 



r(") = In det (i^ + (4"^ + i/i)7o + im) , (39) 



where m is the quark mass and fi is quark chemical potential. The main difference to 
the model discussed in [39] is that we take into account that the gluonic sector is centre 
symmetric and the only centre sector bias arises from the quark determinant. Thus, the 
sum in fl38|l democraticly extends over all centre sectors without any further bias to the 
trivial centre sector. 



4.2 Cold but dense matter 

The calculation of determinant in fl38p can be performed following the techniques developed 
for the Schwinger model. In 4 dimensions the determinant is more severe UV-divergent 
than in 2 dimensions but the effective theory only addresses the low-energy modes below 
a certain physical energy scale. In this approach, the cutoff scale A is finite and acquires 
a physical interpretation. Here, we are only considering temperatures which are small 
compared to this cutoff scale which allows us to effectively set T/A — )• 0. Using a cutoff 
function provided by the Schwinger-proper time approach (see appendix lA.ip . we find: 



r(") = ro(/3,L,A) + ril(/3,L,/x) (40) 
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with cutoff-dependent and cutoff-independent contributions 

ro(/3,L,A) = v^/3A - erf (E(p)/A)] , (41) 

p 

r£i(/3, ^, /i) = 2 5^ { In (1 + e-''(^(^)+'^)) + In (l + < e-''(^(^)"^')) } , (42) 
p 

where we have introduced the 1-particle energy by 

E{p) = v^m^ + p2 ^ p = ?:^(^ni,n2,n^Y' , rik : integer. (43) 

The overall factor 2 in the 4-dimensional case (as compared to the 2d Schwinger model) 
arises from the two spin-orientations of the quarks. We point out that the cutoff dependent 
part Tq is independent from the centre sector number n and the chemical potential while 
those parts T^^J^ which do depend on /i and n are UV finite. 

Let us now consider the case of low temperatures (and n > 0) for which we can neglect 
the anti-quark contributions to T^^J^: 

/3m > 1 Q-P{E{v)+f^) _^ 0. (44) 

For such low temperatures, even the mesonic type excitations can be neglected, and the 
partition function is approximately given by 

^= ,2 



n=l p 

Let us now specialise to = 3, and expand the brackets. We obtain terms such as 



In order to evaluate the sum over the centre sectors, we use 

3 3 ^3 



n=l n=l n=l 



Hence, the centre sector sum eliminates all coloured excitations from the partition function 
by projecting onto states with vanishing iV-ality thus making manifest the confinement of 
colour. For fi < 6m, we obtain the simple result 



3 e^o 



1 + 

Pl,P2,P3 



J2 exp|-/3(E(pi) + E(p2) + ^(p3)-3/i)} 



(46) 



where in the momentum sum only two of the three momentum can be equal. Since our 
model does not include any binding between the quarks, the sum E{pi) + E{p2) + E{ps) 
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Figure 6: Baryon number as a function of the chemical potential for the confining quark 
model (FEC) in comparison to the standard Fermi gas result, left panel. The normalised 
partition function in f l55|) . right panel. 



can be interpreted as the energy of the "baryon" in the context of the present model. For 
chemical potentials smaller but close to the constituent quark mass m, the thermal energy 
density, i.e., 

E..eUT) = , (47) 



behaves like 



-E'therm(T) ^ 3m expj- ^^^^ I , fi<m 



The first excitations in the model are baryonic ones. This result should be compared with 
the case of the "frozen" centre sector (which, in our simple case, is the free quark model): 



El]::miT) ^ m exp{- , /i < m . 



Here, the mass gap is provided by the constituent quark mass, and the first excitations 
which are encountered by increasing the temperature stating at T = are quark excitations. 
As in the Schwinger-model, the centre sector sum solves the Silver-Blaze problem. 



4.3 Condensation of centre dressed quarks 

The thermodynamical quantities of the present quark model can also be calculated exactly 
by evaluating e.g. f H2|) by numerical means. The constituent quark mass m here sets the 
fundamental energy scale. It turns out that the approximations leading to ( l46l) are only 
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valid for rather low temperatures, i.e., /3m ^ 100. In this subsection, we study intermediate 
temperatures and several spatial extensions such as 

/3m = 10 , Lm = 3. . . 10 . (48) 

We are interested in the baryon number Qb which is accessible by 



where the centre sector weights are given by 



(49) 



U^n = '"P^^^"^^, . (50) 

Note that for /i 7^ the weights can be complex and therefore evade a straightforward 
interpretation as probabilistic weight for a given centre sector. An inspection of (H2|) shows 
that 

= u*, u,eR, E^^P^^dci) e M . (51) 

n 

We can always compare our findings with those from a free Fermi gas model for which we 
have 

ui = , U2 = , cus = 1 (Fermi gas) . (52) 

The calculation of the centre sector weights is delicate since the effective actions r„ gener- 
ically are large numbers. To facilitate this calculation, we introduce the subtracted actions 
by 

f „ = r„ - Fmax, Tmax = max (Re F„ j Vn . (53) 
It is easy to check that the centre sector weights are actually independent of this shift: 



E„<=-p{r!,l} E„<=xp{rn} 



Figure El left panel, shows our finding for the Baryon density as a function of the chemical 
potential /i for a small spatial volume: 

mL = 3.5 . 

For fi/m ^ 1.3, we observe that the baryon density strongly peaks for the FEC model. This 
peak is absent in the standard Fermi gas model. The reason for this peak are cancellations 
between the centre sector which nullify the partition function. To see this, we introduce 
the normalised partition function by 



7sub 



exp{rje-r— = Eexp{rSl}. (55) 
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Figure 7: Ratio between the baryon density of the FEC model and that of the standard 
Fermi gas model for two spatial sizes L, left panel. The real part of the centre sector 
weights UJ1...3, right panel. 

Note that Fmax £ ^ implying that a zero of Z^^^ goes along with a zero of the full partition 
function Zq . The quantity Z^^^ also appears in the denominator of the sector weights cuj. 
Figure El right panel, shows the normalised partition function Z^^^ as a function /i. We 
observe for mL = 3.5 that Z^^^ almost vanishes for the chemical potential corresponding to 
the peak position in the baryon density. Increasing the spatial size, we find the near- zero 
regime of Z™*^ is lifted and the peak of the baryon density diminished. For an even number 
of colours, the rise of the baryon number due to a vanishing partition function was called 
Fermi-Einstein condensation (FEC) [321110] • For small enough spatial size, we observe that 
FEC can also occur for an odd number of colours and thus in QCD-like theories. 

Let us now study the baryon density for the bigger systems with mL = 5 and mL = 10. 
Figure [71 left panel, shows the baryon density of the FEC model normalised to that of the 
Fermi gas model. We observe that for small values of fi the baryon density is suppressed 
compared to that from the Fermi gas model. For mL = 10, this suppression ceases for 
fi/m ^ 0.9 and the FEC baryon density equals that of the Fermi gas model. We see here 
confinement at work: while in the Fermi gas model, the baryon density rises with increasing 
H by exciting single quarks into the system. In the FEC model at small values of /i, the 
only way to increase the baryon density is to excite a baryon with mass 3m (in our model). 
For n ^ 0.9m, deconfinement sets in and density rises further on by adding quarks to the 
system. This interpretation is corroborated by an inspection of the centre sector weights cuj 
in figure [71 right panel. Note that corresponds to the trivial centre sector (23 = 1). For 
small values of n, the real part of the weights are roughly the same and equal ~ 1/3. This 
indicates that centre sector transitions frequently occur wiping coloured states from the 
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Lm=10 




Figure 8: Thermal energy E{T) as a function of the inverse temperature /3 for the FEC 
model and the free Fermi gas (left). Ratio between the thermal energy densities of the 
FEC model and the free Fermi gas (right). 

partition function thus installing confinement. For fi > 0.9m, the weights of the non-trivial 
centre sectors, i.e., 001^2, vanish and only the trivial centre sector significantly contributes. 
Hence, the FEC model migrates into the Fermi gas model. 

4.4 Deconfinement at finite temperatures 

Let us now focus on the case of vanishing chemical potential, i.e., /i = 0. Our aim here 
will be to explore the phenomenology of the centre sector transitions far from the dense 
regime. The quantity of main interest is the thermal energy density Ethcrm(T) in (147|) . In 
the FEC quark model, we obtain: 

-E'therm(T) = ^ Un ^ E{p) 
n p 

The centre sector partition functions r^g]^(/3, L, = 0) are real and positive functions of 
the temperature T = 1//3. The centre sector weights Ui, i = 1 . . .3 satisfy 

3 

cji G M, cui > 0, ^ = 1 , (57) 

i=l 

and, thus, can be interpreted as the probability with which each sector contributes to 
e.g. the thermal energy in (1561) . For later reference, we also quote the thermal energy of a 



(56) 
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Figure 9: Probabilistic weights for the realisation of the centre sector n. 
free Fermi gas: 

ElZJT) = 2j2E{p)-^^^. (58) 



p 



Note that antiquarks contribute in the same way as quarks. Hence, the factor 2 in ( 158|) . 
For first insights, we consider the case of small temperatures: 

/3 m < 1 . 

The partition function fl38p can then be approximated by 



n=l p 

3e^» 

p 



= E n + < e-^""^^^] '[l + zn e-''^^^)] ' (59) 



2 



[l+ 2^e-'^2^(*')] . (60) 



For /i = 0, the most important contribution to the low temperature partition function 
arises from mesonic states. Figure El left panel, shows the thermal energy for the FEC 
model in comparison to the free Fermi gas. At large (3 (low temperatures), the FEC thermal 
energy is suppressed like exp{—/3 2m} since only mesonic excitations are possible by virtue 
of confinement. This is in contrast to the free Fermi gas where the contribution from single 
quark states yields a suppression only of order exp{— /3 m}. 
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Dividing the thermal energy by the spatial volume yields the thermal energy density: 

CFEC = Ethcrm/L^, Cfree = -^thcrm/-^'^- 

Figure [HI right panel, shows the ratio of thermal energy densities of the FEC model and 
the free Fermi gas. Around T 0.2, we observe that the FEC energy density rapidly 
approaches that of the Fermi gas. We attribute this behaviour to a deconfinement phase 
transition. Indeed, figure [9] shows the centre sector weights Wj. For small temperatures 
T ^ m, we observe that all weights roughly equal 1/3 indicating an equal contribution 
from the centre sectors to thermodynamical quantities. For T > 0.2m, we find that the 
non-trivial centre sectors cease to contribute implying that the FEC model turns into the 
free Fermi gas theory at high temperatures. 



4.5 The phase diagram from confinement 

We are now considering both finite temperatures as well as a non-vanishing chemical po- 
tential. In order to trace out the phase diagram, we will employ the centre weight U3 G M: 
for a vanishing chemical potential, all weights Ui are real numbers which sum to unity. At 
high temperatures, we have observed in the previous subsection that — 1 while Ui^2 
thus showing confinement. For large non- vanishing chemical potentials, the latter is still 
true: close to 1 signals the transition of the FEC model to the free Fermi gas. Thus, an 
inspection of u-s{fi,T) maps out the deconfinement region in the phase diagram. 

Note that the oo^ < 1 only strictly holds for /i = where the weights Ui enjoy an interpre- 
tation as probabilities. Though in the deconfinement regime will approach 1, there is no 
need for /i 7^ that U3 is bounded from above by 1. In fact, we have seen for small spatial 
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extensions rriL < 2.5 that Fermi-Einstein condensation occurs with potentially large values 
of u-s (see figure [7]). 

Figure fTOl shows our colour coded result for uj3{fi,T). In the large volume limit, e.g. for 
rriL = 15, the result is as expected: there is quite a sharp transition between the hadronic 
regime at low /x and T to the deconfinement regime under extreme conditions. For a smaller 
volume, e.g. mL = 5 (right panel), large values of cus can be observed for low temperatures 
and /i close to the mass threshold. This is the regime for which the partition function 
cancels to a large extent due to centre sector transitions. 

5 Centre sector transitions in the gauged Ising model 



Above using the Schwinger model as well as a SU(3) effective quark theory, we have stressed 
the phenomenological importance of transitions between the centre sectors of the theory 
under investigations. It remains to show that these centre transitions do occur if dynamical 
matter is present: since this matter explicitly breaks centre symmetry, transitions between 
centre sectors could be prohibited in the infinite volume limit. With "matter", we here 
address any dynamical fields (in addition to the Yang-Mills gauge fields) which transform 
under the fundamental representation of the gauge group. An important example is the 
theory of strong interactions, QCD, which is an SU{3) Yang-Mills theory with e.g. three 
(light) fiavours of quark matter. 

In the remainder of the paper, we discuss theories at zero chemical potential and search 
for potential transitions between the centre sectors. Since the Schwinger model does not 
possess genuine phase transitions due to its 2-dimensional nature, we start the consider- 
ations with the ^2 gauge theory coupled to Ising matter in 3-dimensions. This theory is 
easily accessible by means of Monte-Carlo simulations, and offers the possibility of phase 
transitions. We will focus on two competing scenarios which both sketch a different picture 
of the realisation of the centre symmetry: 

(i) The explicit breaking triggers a spontaneous breakdown of centre symmetry. Centre 
sector tunneling does not take place. Observables only receive contributions from the 
trivial centre sector. 

(ii) The explicit breaking of centre symmetry breaking does not bring an end to the 
centre sector transitions. Observables are still averaged over the different sectors. 
There is, however, a bias towards the trivial sector. 

Note that, in scenario (ii), the matter sector is a correction to the gauge sector of pure 
Yang-Mills theory. The Polyakov line expectation value is only non-zero since the lattice 
configurations are biased towards the trivial centre sector. Centre symmetry does also break 
spontaneously at high temperatures though the critical temperature can be quantitatively 
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different from the pure Yang-Mills case. Note that the idea of a spontaneously broken 
symmetry which is also explicitly broken is a useful concept: in QCD, chiral symmetry 
is explicit broken by the current quark masses. Additional spontaneous chiral symme- 
try breaking assigns the pions the role of "almost" Goldstone bosons and explains their 
particular role in the hadron spectrum. 

5.1 Phase diagram of the gauged Ising model 

We are here going to discuss the Z2-gauge theory minimally coupled to Ising spins ax G 
{±1} to ensure Z2 gauge invariance. The model has been considered as the most simple 
gauged Higgs theory, and it has been shown that its phase diagram has the same qualitative 
features than e.g. the SU(2)-Higgs theory [^81 - 152] . As familiar from lattice Yang-Mills 
theories, the gauge fields are represented by the links Z^{x) G {±1}- The partition function 
of the theory is given by 

Z = jvZ^ Va exp{Sz}, (61) 

Sz = A X] Pf,u{x) + K^cr^Z^(x) cr^+^ , (62) 

X,fl>U X II 

where P^y{x) is the usual plaquette 

P^y{x) = Z^{x) Z^{x + /i) Z^{x + u) Z^{x) . (63) 

One easily verifies the invariance of the partition function under gauge transformations 
(fi,. G {±1}): 

a = Vl{x) a^, Z^{x) Zj^(x) = fi(x) Z^{x) fi(x + /i) . (64) 

For K = 0, the Ising matter decouples from the gauge sector, and we are dealing with a 
pure 1j2 gauge theory. It is well known that this theory confines static centre charges, and 
the static potential between two static charges is linearly rising with their distance. The 
Polyakov lines is an order parameter for confinement. For non- vanishing but small values 
K, the dynamical Ising spins can screen the static centre charges, and we obtain string 
breaking: the linearly rising static potential at large distances flattens when the potential 
energy is sufficient to create a dynamical matter pair. Hence, the phase diagram as a 
function of k and /3i is expected to show the same features as that of the SU(2) Higgs 
theory. 

Figure [HI left panel, summarises this diagram: due to the lack of a local order parameter, 
a cross-over region from the string breaking regime to the so-called Higgs phase is expected. 
We verified these expectations by calculating Polyakov line expectation values using a 20^ 
lattice. Both, the update of the links Z^{x) and the matter fields ax have been done by 
standard heat bath techniques. We calculated the Polyakov expectation value for 204 x 201 
different values in the [5i-k plane. Our numerical findings are colour-coded and summarised 
in figure [TTl right panel. Both, the cross-over regime as well as the finite size transition is 
clearly visible. 
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Figure 11: Sketch of the phase diagram (left); Polyakov line expectation value of the gauged 
Ising model as function of /3i and k. 



5.2 The interface tension of the trivial centre sector 

Let us briefly discuss the "empty" vacuum of the gauged Ising model. The discussion in the 
context of pure Yang-Mills theory in subsection 12.11 can be straightforwardly transferred to 
the case of the 7^2 gauge theory: The gauge inequivalent empty vacua are characterised by 
the values of the homogeneous Polyakov lines in each direction. Thus each of the Polyakov 
lines Pi, P2 and P3 takes values ±1, there are 2^ = 8 states which we need to consider. 
The ^2 centre transformation 



Zjx) (-DZjx) 



V X., yu 7^ z/, 



flxed 



(65) 



is discrete, changes the sign of the Polyakov line and therefore maps one empty vacuum 
state to another. 

Before we proceed to consider the centre sector transitions within the gauged Ising model, 
we here discuss centre interfaces in the particular sector with all Polyakov lines trivial. For 
this "empty" vacuum state a gauge can be found where all links are one: 



iPuP2,P3 



:i,i,ii 



Zf^{x) = 1 V X, /i . 



(66) 



The gauged Ising model then collapses to the standard Ising model with ferromagnetic 
bonds only. 

If is particular 2d array of spin at xq = 0, i.e., the so-called m-state and if is 
a centre copy of this state for which all spin are reflected, we would like to investigate 
the overlap of the true vacuum state |V^o) with either Itp) and If centre symmetry is 
realised in the Wigner mode, the vacuum is centre symmetric yielding an overlap in both 
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Figure 12: Partition function with anti-periodic boundary conditions rewritten as a parti- 
tion function with a centre twist. 

cases. If the centre symmetry is spontaneously broken, the overlap vanishes for one of the 
matrix elements. To investigate the spontaneous breaking of centre symmetry, we thus 
study: 

^ (>|exp{-ff/r}|^) 
^ (^|exp{-i//T}|^> ' ^'^^ 

where H is the Hamilton operator and T is the temperature. For sufficiently small tem- 
peratures, the exponentials in the latter equation project onto the ground state I'lpo), and 
we obtain: 

^ (V|^o) (V^o|V^) _ f 1 Wigner-Weyl realisation 

I (■?/'o|'?/') 1^ I spontaneous symmetry breaking. 

In the functional integral approach, x is given by the ratio of two partition function: the 
partition function with the spins cTx obeying anti-periodic boundary conditions constitutes 
the numerator in (!67|) . while the standard partition function with period spins is in the 
denominator. A connection to the centre symmetry can be established by "pushing out" 
the centre element to the links: 

U^{x) ax+^ = (-1) U^{x) ax+^ . 

The net effect is illustrated in figure [T2l 

7 7 

anti— periodic -^twist I R(X\ 

X ~ y ~ ~^ • V^'^J 
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We expect that the transition between centre sectors are exponentially suppressed with a 
slope given by the size of the centre interface. For a A^^ x A^^ lattice, the minimal surface 
of the centre interface is iV^. It is therefore convenient to introduce the interface tension 
a by 



For small values k, x can be directly calculated by the so-called strong coupling expansion, 
i.e., the Taylor expansion with respect to n. For large values of k, x can be obtained by 
means of a duality transformation. In 3d Ising model is dual to the 3d gauge theory, 
and X appears to be the expectation value of the maximal spatial Wilson loop (details 
will be presented elsewhere). Altogether, we find: 



For a fixed aspect ratio Nt/N and in the thermodynamic limit A^ — 00, we find that 
the theory is in the disordered phase for small i.e., x implying that centre sector 
transitions do occur frequently. For large k, the interface tension is independent of Nt and 
large showing that spontaneous centre symmetry breaking occurs in this case. 

For intermediate values of k the interface tension must be obtained by numerical means. 
The ratio of partition functions, i.e., Xi can be numerically calculated in an efficient way 
using the so-called snake-algorithm [53]. Our numerical results are shown in figure [131 In 
the ordered phase at k = 0.20, we observe an exponential decrease of the interface tension 
cr with A^ in line with result from the leading order Taylor expansion in k. For the ordered 
phase at K = 0.25, we confirm that the interface tension is indeed independent of the 
volume. 

5.3 Centre sectors and parametric transition 

While in the previous subsection, the link variables have been frozen to the trivial centre 
sector, we now consider the transition element x of the gauged Ising model where all links 
Z^{x) are dynamical. Defining the twisted link variables by 




(70) 




for K , 
for K > 1 . 



(71) 




Vxi, X2 ; xs fixed , 
else. 



(72) 



the partition function of the twisted partition function is given by 




(73) 



(74) 



Given the invariance of the measure and of the plaquette. 
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Figure 13: The interface tension of the 3d Ising model as a function of the system size N 
for interaction strength k = 0.20 (disordered phase) and k = 0.25 (ordered phase). 



one easily obtains 

-^twist = -Aperiodic ^ X = 1 • (75) 

In the gauged Ising model (as well as in the Schwinger model and QCD), the centre twist 
is part of the gauge field average. It is therefore, strictly speaking, impossible to consider 
a spontaneous breakdown of the centre symmetry. Note, however, that the matter fields 
break the centre symmetry explicitly. Depending on the coupling strength k and the 
temperature, Nt respectively, this explicit breaking can be small. Actually for small k, we 
are going to show that at low temperatures the explicit breaking is rather irrelevant while 
above a critical value for the temperature, the explicit breaking is strong in the sense that 
it does not makes sense to even consider an approximate centre symmetry. We call this a 
parametric transition. 

Let us firstly illustrate the explicit breaking in the gauged Ising model. To this aim, 
we integrate over the Ising spin for a given link distribution Z^{x). The resulting 
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Figure 14: Polyakov line expectation values of the gauged A^^ Ising model as a function of 
the system size (left). Same quantity for a 6 x A^^ lattice as a function of /3i (right). 

probabilistic factor contains contractible Wilson loops W[Z] as well as the Polyakov line 
P[Z] (for simplicity we here consider infinite space) 

I Va eMSz\Z]} = /^^^''^ Wc[Z] + P[Z] + O . (76) 

C Xi,X2 

The sum over the closed loops C is centre symmetric while the terms containing the 
Polyakov loop explicitly break the symmetry. 

In following, we study the theory by means of Monte-Carlo methods. We used a heat-bath 
update for the link and matter fields. Since the matter action describing the gauge field 
matter interactions is local, we used a generalised Luescher-Weisz algorithm for the 

calculation of the Polyakov line expectation value. 

We are now considering moderate values for k, let us say n < 0.25. We have already 
obtained (see figure [TT]) that for small Pi the Polyakov line is small. Increasing Pi above 
a critical value, the Polyakov line expectation value (P) rapidly takes large values. In 
fact, we find that, below the critical value, (P) is strictly vanishing in the infinite volume 
zero temperature limit (for fixed and finite aspect ratio) although (P) 7^ for any finite 
system. This is illustrated in figure [TU left panel: the graph shows (P) for a lattice as 
a function of A^. At least in the string-breaking phase for sufficiently small k, we observe 
an exponential decrease of (P) with increasing A^. In the Higgs phase (here for Pi = 0.6 
and K = 0.30), we still observe a slight decrease of (P) with the system size. Whether the 
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Polyakov line expectation finally vanishes for large cannot be inferred from the present 
data. Further investigations are left to future studies. 

To monitor the deconfinement cross-over at finite temperatures, we have calculated (P) 
for a 6 X A^^ lattice as a function of /3i for the aspect ratios 4, 5, 6, 7. Our numerical result 
is shown in figure [TU right panel. We observe a maximum of (P) for /3i = /3c ~ 0.75(5) 
and very little dependence on the system size N. We argue that Pi < Pc, the explicit 
breaking of centre symmetry is weak, the physics of the regime might be similar to that 
of the symmetric theory for k = 0, and that perturbation theory with respect to k, yields 
reliable results. On the hand for Pi > P^, centre symmetry breaking is strong and any 
similarities with the symmetric theory are lost although we always have x = 1 in f TTSj) . 
This is the parametric transition. 

6 Centre-sector transitions in the SU (2) Higgs theory 



6.1 String breaking 

In order to study centre sector transitions in a more realistic, i.e., more QCD-like, the- 
ory, we will study a theory where the" empty vacuum" features flat directions, where this 
vacuum symmetry collapses to the discrete centre symmetry upon the inclusion of quan- 
tum fluctuation and where the trivial centre-sector is favoured through the matter sector. 
In addition and for the first time in this paper, we will consider 4 Euclidean space-time 
dimensions. 

The most simple such theory is a SU (2) gauge theory with a scalar field (say Higgs) in the 
fundamental representation. This theory is accessible by means of lattice gauge theories 
with large statistical accuracy. The partition function is given by 

Z = jvUV<p P0t exp{5wii + 5Higgs} , (77) 
^Higgs = Re0^(x)f/^(x)0(x + ^) - \}^4,\x)<P{x) + \[ci>\x)<P{x)f . 

The parameter k quantifies the interaction strengths of the scalar fields with the gauge 
sector, while A gives rise to a quartic Higgs self-interaction. 

The SU(2)-Higgs theory shares the so-called string breaking with QCD: at low tempera- 
tures and sufficiently weak gauge-matter interactions, the static quark anti-quark potential 
linearly rises at intermediate distances between the static quark antiquark pair while it fiat- 
tens at asymptotic distances. A popular picture assumes that a colour electric flux tube 
still forms between a static quark antiquark pair and that this string breaks due to the 
creation of a Higgs anti-Higgs pair. Once the string is broken, there is little penalty in en- 
ergy to increase any further the distance between the (screened) heavy quarks. The static 
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Figure 15: For the SU(2) Higgs theory (A = 0.1) the static quark antiquark potential for 
a 20^ lattice, /3wii = 2.2 and several values for n, 8160 configurations per set. 
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potential approaches a constant. This also implies that the Polyakov line expectation is 
non-zero (even at low temperatures) and cannot serve any more as an order parameter. 

The phase diagram has been qualitatively described in figure [HI On a more quantitative 
note, we have calculated the static potential by means of lattice gauge simulations. In 
order to ensure good ergodicity, we used the local hybrid Monte- Carlo update for the link 
fields as well as the scalar fields. In this paper, we work with A = 0.1 and /3wii ^ [2.2 . . . 2.5] 
which roughly corresponds to the scaling window, at least, at small to moderate values k. 

To calculate the static quark potential, we here followed the two channel approach devel- 
oped in [56],[57]. Using a 20"^ lattice, our results for the Wilson coupling /Swii = 2.2, fixed 
quartic coupling A = 0.1 and several values of k is shown in figure [T51 For each potential, 
8160 independent lattice configurations contributed the relevant expectation values. For 
small values of k, e.g. k, < 0.336, the string breaking scale can be intuitively set by an 
inspection of the graph: for k = 0.31, we would set r;, ~ 7 a. For larger values of k, a grad- 
ual transition to the Higgs phase sets in, and a string breaking scale is hardly defined. It 
seems, however, that this scale cannot be pushed to small values by just altering the Higgs 
hopping parameter k. If we fit the static potential for r < to V{r) = oq + CTint?' + cli/^, 
we define the so-called intermediate string tension cxint. For k, = 0.31, our estimate is 



K values close to the transition line, e.g., for k = 0.336, the string breaking scale is difficult 
to define. Here, the picture of a well defined short string might loose its validity. 

6.2 Volume dependence of the Polyakov line expectation value 

Integrating over the Higgs fields for a fixed link background yields an effective action for 
the links which is not anymore centre-symmetric. In order to get a first impression of 
the amount of explicit breaking we consider the normalised probability distribution for the 
spatial average of the Polyakov line. 



where the expectation value is with respect to the partition function ( fT2l) and P[U] is the 
spatial average of the trace of the Polyakov line: 



Since in the low temperature phase the probability distribution of Polyakov line V is 
basically given by the Haar measure of the gauge group, the above distribution is divided 
by the reduced Haar measure distribution. 
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Figure 16: Ratio R{p) of the Polyakov line distribution function over the Haar measure 
distribution for three different temperatures in pure Yang-Mills theory (left) and for the 
SU(2) Higgs theory on the 20^ x lattice (right). 



to extract any effects beyond the trivial distribution. 

Let us first study pure Yang-Mills theory without Higgs matter. Our result is shown in 
figure fT6l left panel. At zero temperatures, i.e., Ng = Nt = 20, for /3wii = 2.5, this ratio 
is basically independent of p. We increase the temperature T = 1/Ntahj decreasing the 
extent of the lattice in time direction. For A'^j = 6 at /3wii = 2.5, pure Yang-Mills theory 
is still in the confinement phase. Here we observe that deviations from the Haar measure 
distribution is still marginal. For Nt = 4, Yang-Mills theory is in the high temperature 
deconfinement phase, and significant deviations from the Haar measure distributions are 
clearly visible for large values for p. Here, a bias towards the centre elements, 1 and — 1, 
is seen. In the infinite volume limit, tunneling between the centre sectors ceases to take 
place leading to a spontaneous breakdown of centre symmetry. It is, however, important 
to notice that the Aj = 4 distribution W{p) is to a good extent symmetric under the 
reflection p — )■ —p by virtue of the centre invariance of the functional integral. This signals 
good ergodicity of the algorithm even in the high temperature phase. 

Let us now consider the SU(2) Higgs theory. We here used a 20^ x Nt lattice and /3wii = 2.5, 
A = 0.1 and k = 0.31 which is somewhat below the critical value for the transition from 
the confining phase to the Higgs phase. Figure Uni right panel, summarises our findings. 
At very low temperatures, i.e., for Nt = 20, there is hardly any deviation from the Haar 
measure distribution visible. Still in the string-breaking phase at At = 6, there is a bias 
towards positive values of the Polyakov line noticeable. At high temperatures, such as 
for At = 4, above the deconfinement transition, this tendency is strongly amplified: the 
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Figure 17: For a fixed aspect ratio N/Nt, the dependence of tlie Polyakov line expectation 
on Wilson /3vvii in the SU(2)-Higgs theory k = 0.31, A = 0.1, 2040 configurations per data 
point, (left panel). Same quantity as a function of InT in the Schwinger model (right 
panel) . 



centre sector 2:1 = —1 is even more suppressed while the probability for positive values of 
p reaches large values. In all cases (though less visible for A^^ = 20), the centre symmetry 
U — )■ If^, f l8|9|) is explicitly broken, and the Polyakov line probability distribution is no 
longer reflection symmetric, W{p) 7^ W{—p). 

Apparently, the explicit centre symmetry breaking is much stronger at high temperatures 
than in vacuum. Does this mean that centre symmetry breaking is spontaneously broken 
at high temperatures? To answer this question, we have studied the the Polyakov line 
expectation value (P) for a fixed and finite aspect ratio, i.e., N/Nt, as a function of the 
system size. This can be easily achieved by using of fixed number of lattice points, 16^ x 6 
and 24^ x 6, respectively, and to vary the Wilson /3wii- The logarithm of the temperature T 
in relation to the fundamental renormalisation group scale such as the intermediate string 
tension a is then roughly given by 

In(^) = In ( J" — -) ^ 7i/5wii + constant, 

where 71 is 1-loop Gell-Mann Low coefficient. For a pure SU(2) gauge theory, this coefficient 
would be 7i = 37r^/ll. Figure [T71 left panel, shows (P) as a function of Wilson using 
the Liischer-Weisz method [511IS5] which is easily generalised to include the Higgs field. 
For small values of /3wih i-e., /3wii < 2.4, we roughly observe an exponential decrease of (P) 
with decreasing /3wii- At high values, i.e., /3wii ~ 2.4, (P) increases with increasing at a 
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modest pace. Whether in this regime the dependence of (P) is still exponential (just with 
a much smaller slope) or whether the characteristic of this dependence has fundamentally 
changed cannot be concluded given the present set of data. 

In pure SU(2) Yang-Mills theory and a 16"^ x 4, we would expect the deconfinement phase 
for /Swii ~ 2.4. It is tempting to conclude that the change of the dependence of (P) on 
in this regime signals a spontaneous breaking of centre symmetry. We here argue that this 
conclusion is premature: figure ITTl right panel, shows (P) as a function of In T/m^ (where 
is the induced photon mass which sets the scale) for the same fixed aspect ratio. Both 
curves are strikingly similar. While there is certainly no spontaneous breakdown of centre 
symmetry in the Schwinger model (this would lead to a Silver-Blaze problem as we pointed 
out in section |3]), the large values of (P) merely indicates an overlap problem if we wish to 
address this model by means of Monte-Carlo methods. 

6.3 Order parameter for centre-sector transitions 

If we consider for the moment a given lattice configuration {U} in pure Yang-Mills and 
its centre copy {^U}, then these configurations are degenerate in action. If {U} is an 
"empty vacuum" state, i.e., all contractible Wilson loops on the lattice yield the unit 
element, a path in configuration space can be found which deforms {U} into {^U} without 
changing the action. For a generic lattice configurations, this no longer the case, and the 
crucial question is whether transitions between centre-sectors occurs at all in the infinite 
volume limit. The situation is aggravated if dynamical matter (transforming under the 
fundamental representation of the gauge group) is present which bias the trivial centre 
sector. We have partially answered this question by studying the volume dependence of 
the Polyakov line: in the string-breaking phase, the bias becomes negligible in the infinite 
volume limit, centre transitions do occur and the Polyakov line averages to zero; in the 
Higgs-phase explicit centre-symmetry breaking is strong and independent of the volume. 

Our aim here will be to construct a sort of order parameter which is sensitive to the centre 
transitions. We do not expect that such an observable is built from local field operators 
and that it is an order parameter in the strict thermodynamics sense. It must necessarily 
be a non-local object which, however, nevertheless signals whether swapping the centre 
sectors can occur. To this aim, we divide the 3-volume into two parts of equal size, Vl and 
Vr and define the spatial average of the Polyakov line over each of the volumes: 



It is straightforward to assign a centre sector to each of the sublattices by the mapping 





C(P) = n. 



n : arg(P) 
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arg(P)=^6]0,27r] 
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Figure 18: The tunneling coefficient for pure Yang-Mills theory and the SU(2) Higgs theory. 

Let us then constrain the configurations in such a way that the left hand part of the 
universe belongs to centre-sector n while the right hand part has a reference to centre 
sector m. The ratio of partition functions between the constrained and the un-constrained 
theory defines the free energy Fnm for the mixing of the centre-sectors: 

:= exp{-F„^/T} = ^ j VU Vct>Vct>U(n,C{PL)) 5(m,C(P^))e^ (84) 
N = j VUVcj) Vct)^ e^ . 

The matrix pnm can be interpreted as the probability to find sector n in the left half and 
sector m in the right half of the spatial universe. If the centre-symmetry is weakly broken 
in the string- breaking phase (or not broken at all in the infinite volume limit), the free 
energy Fnmi n ^ m is finite (or might even tend to zero for an increasing system size), 
while the free energy diverges if the centre symmetry is strongly broken broken as e.g. in 
the Higgs phase or at high temperatures. This is basically due to the fact that the whole 
universe belongs to one centre sector, and configurations with n ^ m have a very low 
probability. 

Still measuring the free energy Fnm is hardly an easy task. Alternatively, we consider 
another intuitive measure for centre sector transitions: the probability r that the volumes 
Vl and Vr belong to different sectors. We call r the transition coefficient. It is directly 
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related to the free energy by 



r = ^Pnm = 1 - ^Vnn = 1 - ^ exp{-F„„/T} . (85) 

n^m n n 

In the infinite volume limit in the string-breaking phase, explicit centre symmetry breaking 
can be neglected implying 

_ 1 _ 1 

In the Higgs or the high temperature phase, the whole lattice belongs to one centre sector, 
and we find: ^ 

Pnm ^nm • 

We have numerically estimated the transition coefficient r using a 24^ x 6 lattice, k = 0.31 
and A = 0.1. The parameter setting is such that the theory is in the string-breaking phase 
for a 24^ lattice size. Figure [18] shows our findings for r as a function /3. We do find that for 
/3wii ~ 2.35 centre-sector transitions do occur with high probability while one centre-sector 
is observed throughout the lattice universe for high temperatures, i.e., /3wii > 2.35. 



7 Conclusions 

For SU(N) Yang-Mills in 4 dimensions on a torus, we started to investigate the "empty 
vacuum", i.e., the set of configurations for which all holonomies calculated for any con- 
tractible loop yields the unit element of the group. All these configurations produce zero 
field strength everywhere. We found a continuous set of gauge in- equivalent configurations 
related by a symmetry transformation. Including quantum fiuctuations, the degeneracy 
of these states is lifted. The symmetry collapses to the well-known Z^r centre symme- 
try which divides the gluonic configurations into centre sectors. Transitions between the 
centre sectors turned out to be the key ingredient for confinement Including dynamical 
matter which transforms under the fundamental representation of the group, this centre 
symmetry is explicitly broken. Our central working hypothesis to start with was that tran- 
sitions between centre sectors still take place in the so-called hadronic phase and that 
these transitions only cease to exist at high temperatures when the centre symmetry is 
also spontaneously broken (on top of the explicit breaking). 

Before corroborating this picture in the sections |5] and |6l we studied the phenomenological 
impact in the Schwinger-model, since it allows for explicit analytical solutions, and in an 
SU(3) quark model for the sake of its relevance for QCD. For an even number of colours, 
it was firstly pointed out in [391110] that quarks acquire periodic boundary conditions 
in some of the centre sectors. For a finite chemical potential, the quarks then might 

■^although they do not explain the confinement energy scale of several hundred MeVs for QCD 
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undergo condensation due to Bose statistics. In analogy, this has been called Fermi Einstein 
condensation (FEC). We traced out the roots of FEC in the Schwinger model on the 
torus as finite chemical potential. We found that the phenomenological importance of the 
centre transitions is the solution of the silver blaze problem: centre transitions wipe any 
dependence of the partition function on the chemical potential as it must be since the 
physical states of the model carry no net baryon number. 

We then studied an effective SU(3) quark model with constituent quark mass m which only 
interactions are with the gluonic background field specifying the centre sector. We verified 
that this very simple model already confines quarks: e.g., we considered the thermal energy 
density as a function of (low) temperatures and found that its lowest excitations are mesons 
of mass ~ 2m rather than quarks with mass ~ m. By studying the centre sector weights 
(provided by the model), we found that the centre sectors democraticly contribute to the 
partition function in the hadronic phase while under extreme conditions basically only the 
trivial centre sector contributes implying that the model merges with the Fermi-gas model 
in the quark gluon plasma phase. Using the centre weight of the trivial centre sector as an 
order parameter, we were able to map out the phase diagram of the model as a function 
of the chemical potential and the temperature. We found that FEC is important at low 
temperatures and intermediate values of the chemical potential. 

For FEC to happen in dense QCD, the central question is whether centre sector transitions 
do take place despite of the explicit breaking by matter fields. This question can be studied 
for vanishing chemical potential where Monte-Carlo simulations are readily available. Using 
a 3-dimensional gauge theory with Ising matter, we firstly studied the tension of a centre 
interface when the model is forced into the trivial centre sector (i.e., the sector for which the 
model is identical to the standard ferromagnetic Ising model). While the interface tension 
vanishes in the infinite- volume zero-temperature limit, it is finite at high temperatures 
when centre symmetry is spontaneously broken. Giving up the artificial constraint which 
ties the gauge theory with matter to the standard Ising model, we find that interface 
tension vanishes for any finite volume. Nevertheless, a detailed study of the Polyakov 
line expectation value reveals at rather small value (though non-zero) in the would-be 
confinement phase, and large values in the Higgs phase. 

A theory which is more relevant for QCD, is the SU(2) gauge theory with dynamical 
matter. We have chosen a scalar Higgs field since it allows large statistic Monte-Carlo 
simulations (facilitated by e.g. the generalised Luescher-Weisz method) and at the same 
it explicitly breaks centre symmetry as quarks do. We firstly identified the region of 
the coupling space where the theory realises string breaking well within the size of our 
lattice. We then studied the extent of explicit centre symmetry breaking by calculating 
the Polyakov line distribution function. At high temperatures, the explicit breaking is 
largely amplified by an additional spontaneous breakdown. A study of the temperature 
dependence of the Polyakov line showed quite a similar behaviour as in the Schwinger 
model. This renders cumbersome any conclusion on the spontaneous centre symmetry 
breakdown since it is clearly absent in the Schwinger model. This calls for a new type of 
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order parameter which is designed to be a Litmus paper for centre transitions. For such an 
order parameter, we here proposed to map the Polyakov hne to the centre sector for each 
half of the spatial lattice universe separately, and ask for the probability that a particular 
lattice configuration belongs to different centre sectors in each half. If this is the case, 
the theory certainly undergoes transitions between the centre sectors. Using this order 
parameter, we find clear numerical evidence in the SU(2) Higgs theory for a "hadronic" 
phase at small temperatures and a de-confined phase at high temperatures. 

In conclusion, we found evidence that dynamical matter does not prevent centre sector 
transitions in QCD-like theories in the low temperature "hadronic" phase. These tran- 
sitions only cease to exist at high temperatures when centre symmetry is also broken 
spontaneously. We highlighted the phenomenological impact of these centre transitions: in 
the Schwinger model, they solve the Silver-Blaze problem, and in an SU(3) quark model 
they lead to a new phase featuring Fermi Einstein condensation (FEC) for cold but dense 
matter. The FEC mechanism only uses fairly robust assumptions on the realisation of 
centre symmetry. Hence, FEC might also be at work in QCD at low temperatures and 
intermediate values of the chemical potential for which quarks are not yet liberated. 
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A The Schwinger-model calculations 

In this appendix we extend known results about the Schwinger model at finite temperature 
T = [29] to the case with chemical potential /i 7^ 0. For related results on the Thirring 
model see [32] • We enclose the system in an interval with length L and impose periodic 
boundary conditions in the spatial direction. 

A.l Schwinger-proper time regularisation 

Any gauge potential with vanishing instanton number on the two-dimensional torus can 
be decomposed as 

27r 27r 
Ao = -di(j) + doX + jho , A, = do(j) + d,X + —h (86) 

with constant toron fields Hq and hi. The determinant of the Dirac operator for massless 
fermions does not depend on A and its 0-dependence follows from the axial anomaly |29] . 
Here we calculate the determinant for arbitrary toron fields /i^ and for a chemical potential 
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/i 7^ 0. For constant fields the Dirac operator possesses plane waves as eigenf unctions and 
its determinant can be written as an infinite product over all admitted momenta, 

det(i4j= n Df,,^im,En) , (87) 

(m.,n)eZ2 

The eigenfunctions are anti-periodic in time and periodic in space such that 

DHAm,En)=(^j^ {m + \--fY + El ^ = ho + '^fi, (88) 

where LE^ = 27r|n — hi\. We introduce the fermionic effective action Th,^ by 

r^,^ = lndet(i4,J (89) 
and obtain in Schwinger proper time regularisation 

poo J 

rA(/x, h)= - V e-^^''-('"'^") ^ V T,UE^, h, /x) . (90) 

/i/A2 S 

I in,n n 

The terms in the last sum are given by the integral 

poo 

rfcr(^,/i,/i) = 2 / deeW{e), (91) 
Je 

where the integrand contains the function 

POO 

WiE) = dsj^ exp{-sD,.^(m,E)} .. (92) 

Now we can do the proper time integration. Since the resulting sum over m is convergent 
we may set T/A = such that 

iy(i?)=e-^'/^' J]/(m,E) with f{m,E) = —-l—-. (93) 

The sum over m can be rewritten with the help of the Poisson resummation formula 
fim, E) = J2 fik, E) , f{k, E)= fdx e^-^^ /(x, E) . 

n k 

After a shift of the integration variable the last integral takes the form 

o2 2Trikho > ■ ^ 



m E) = J . (94) 
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The integral over x is known in closed form and the resulting sum over k can easily be 
computed and yields 



1 d 
2EdE 



{f3E - log (1 + e^-^-^^) - log (1 + } . (95) 



Multiplied with exp(— E'^/A^) this becomes the function W{E) in fPT]) . The /i-dependent 
part is UV-finite and hence we can safely remove the cutoff there leaving us with 

Ta{(3, ft, h) = (3Ea{L, h) + ri(/3, L, /i, h) , (96) 

with a divergent zero-point energy and a finite temperature correction 

ri(/3, L, /i, /i) = E In I (1 + e'^^'^o e-^iE„+^)) + In (l + e"'"^'^" g-'^^^"-^)) } . (97) 

In the zero-temperature limit we must recover the well-known Casimir energy of fermions 
on a circle with circumference L 



(9J 



The Casimir energy is periodic in hi with period 1 and in the last formula we must assume 
hi G [0, 1]. We conclude that the renormalised effective action is 



r(/3, L, /i, h) = -f3EcUL, hi) + Ti{f3, L, /i, h) . 



(99) 



Note that the energies En are proportional to 1/L such that Fi depends only via the 
dimensionless parameter r = on the size L of the system. 

A. 2 6'-function representation and integration over toron fields 

The effective action T = logdet(i(^) can be expressed in terms of the 9- and ?7-function [32] 



1 



det(i4,M) = ^7^0 



rfiiT) 



hi 



1 



(0,ir)e 



hi 



-7 



(0, it) 



(100) 



where 7 is defined in f l88|) and we used the Dedekind eta-function 

r/(ir)=e— /^^nil-^"'^™)' ^=7' 
and the theta function [60] 



n>0 



(101) 



e 



(0,ir) = 5^e- 



■7rr(n+Q!)^+27ri(n+a)7 



(102) 
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which has the product expansion 



' e 



r]{iT 



(^Q j^^ _ g27ria:7g-7rr«2+7rr/12 
oo 



n=l 



Note that the determinant is invariant under large gauge transformations ^ + 1. 

Using the product expansion for the ^-function it follows at once that the effective action 
T{f3, L, fi, h) = \ogdet{i0h,ij.) has the series expansion (!99|) . In the derivation we assumed 
that hi takes its values in the unit interval and this is why the /ii -periodicity of the 
determinant is not manifest. In the thermodynamic limit we find the following expression 
for the free energy density / in Z = exp(— /3L/): 

Jo 

= At [dilog (1 + 6^"''^°-^'^) + dilog (1 + e-■^^^^ho+|3^^\l _ Q04) 

L V / \ /J 

For L//3 — )■ cxo the free energy density does not depend on the constant gauge field hi, 
as expected. When one varies the constant gauge field ho (or equivalently the boundary 
condition in the temporal direction) then one smoothly interpolates between free fermions 
and free bosons or between a Fermi-Dirac and a Bose-Einstein distribution. 
In a gauge theory we must average over all gauge fields and in particular we must integrate 
the fermionic determinant over the constant gauge fields as well. The integral of f llOOp over 
hi can be done explicitly. If m, n are the summation indices in the double sum f llOOl) we 
change summation indices according to m = p + q and n = q. The sum over q together 
with the /ii-integral over [0, 1] turns into an integral over [—00, 00] and yields the simple 
result ( 127|) . The final integration over ho is easily performed and yields the /^-independent 
result (EH]) for the averaged determinant J dhodhi(\et{i(j)f^^^. 

A. 3 Polyakov loops 

Let us finally calculate the expectation values of product of Polyakov loops 

P^{u) = g«'?/'^^°eAo(a;'',u) _ ^2mqho ^-ieq J du<j>{x° ,u)dx'^ (105) 

corresponding to static charges g G Z. In particular P^q = Pg. First we do the integration 
over the harmonics. With the help of f l27|) we obtain 



J dhodhi det (14,^) n ^3'' 



N ^ ^ ^-KTQ^/2+l3nQ^-i(cfxj[x)4>{x) 
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with total charge Q = "^Qi and source 

j{x) = eJ2QiS'i^' -Ui). (106) 

i 

The functional integral over has been done previously in [22] and yields 

'n^^-("^)) = ll{P,^)e-^^'^^''^^^''^\-'-'^^\^ (107) 



with expectation values (133!) for the individual Polyakov loops and periodic potential 

vrm^ cosh {^{L — 2\u\)) 



V{\u\) 



4 



sinh 



m-/L 



\u\< L . 



(lOJ 



The potential energy decreases exponentially fast for large separations of the charges 



V{\u\) 



vrm. 



7 -m^L/2 



(109) 



B Local-hybrid Monte Carlo for the Higgs sector 

For the update the scalar field in accordance to the functional integral (177|1 . we here discuss 
the Local Hybrid-MC (LHMC) scheme. Assume that (pix) is chosen for the update. The 
LHMC Hamiltonian is given by 

1 



H 



(110) 



To set up the LHMC scheme, it is convenient to work with real variables, i.e., vr = R + iS, 
(p = r + is, s, S,r, R E M.. The Hamilton- Jacobi equations of motion are given by 



R 



dH 

'or 

dr 



R. 



. _ dH 

' ~ 'dS 
dSL dSL 



d(j) dcf)'^ 



S, 

s 



ds 



d(f) 90t 



;iii) 

;ii2) 



It is easy to check that the above equations imply H = 0. Combining both equations in 
(Iml) yields: 

TT = R + iS = -2—4, = r + is = vr. (113) 

a0T 

For the action in (177|) . the terms of S'niggs which depend on 0(x) or 0"''(x) give rise to 
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Hence, defining 

B{x) = J2[u,{x)<p{x + ix) + Ulix-^^)<P{x-^l)] , (114) 

we finally find: 

= TT, TT = kB{x) - [l +4A0t(x)0(x)j . (115) 
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